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ABSTRACT 


v 


Analytical  expressions  have  been  developed  for  determining 
coefficients  for  cones  and  spheres.  .  T-he  ranges  Of ~vaiiab.es 
considered  are: 


Cone  half  angle 

Cone  length 

Cone  bluntness  ratio 

Sphere  diameter 

Wall  temperature 

Free-stream  Mach  number 

Altitude  or  approximate 

Free-stream  Reynolds  number 

Angle  of  attack 


8  to  1  5  degrees 
1  to  12  feet 
0  to  0, 3 

3. 75  to  !  2  inches 
1000  to  3000°  Rankine 
5  to  25 

0  to  200, 000  feet 
104  to  1G10 
0  to  20° 


^An  error  analysis  is  presented  in  which  the  resulting  drag 
equations  are  compared  to  wind  tunnel  and  ballistic  range  data, 
and  some  limited  flight  test  data.  'v 


(>Also  presented  are  the  results  of  an  investigation  aimed  at 
determining  the  altitude  range  where  a  conical  or  spherical  re-enir> 
vehicle  will  undergo  boundary  layer  transition. 


This  report  contains  essentially  the  same  technical  information 
as  RAD  TR  65-16  (SECRET)  which  was  prepared  by  the  same  authors. 
However,  the  classified  portions  of  this  TR  have  been  deleted  in 
order  to  permit  a  more  general  distribution  of  useful  data. 
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NOMENCLATURE 


C  Form  at  Chapman -Rube sin  viscosity  coefficient  (also  pressure  coefficiei 

Cp  Specific  heat  of  constant  pressure 

C{  Local  skin-friction  coefficient  based  on  free  stream  conditions 

CDf  Skin -friction  drag  coefficient  (no  interaction) 

CD  Base  pressure  drag  coefficient 

3 

Cp  Forebody,  pressure  drag  coefficient  (no  interaction) 

P 

Cp.  Induced  drag  coefficient 

Wall  temperature  parameter 
Fj  (K)  Induced  pressure  gradient  parameter 
h  Enthalpy 

HT  Total  heat  absorbed  by  body 

X,  K',K"  Induced  drag  functions 
L  Slant  length  of  sharp  cone 

m  Mass  injection  rate 

M  Mach  number 

Pr  Prandtl  number 

p  Pressure 

<;*  Effective  heat  of  ablation 

,,  Dynamic  pressure 

r  Body  radius  of  curvature  (also  recovery  factor) 

Rb  Base  radius 

Nose  radius 
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R„  Reynolds  number  based  on  free  stream  conditions 

X 

Rcx  Reynolds  number  based  on  local  conditions  at  boundary  layer  edge 

T  Temperature 

u  Velocity  in  streamwise  direction 

v  Velocity  normal  to  the  wall 

x  Coordinate  along  sharp  cone  surface 

5*  Boundary  layer  displacement  thickness 

X  Viscous  interaction  parameter 

A  Denotes  induced  values 

a  Angle  of  attack 

rj  Transpiration  factor 

6  Cone  angle  {aemi-vertex  angle) 

0s  Shock  angle 

y  Ratio  of  specific  heats 

p  Density 

p  Viscosity  coefficient 

&  Meridian  Angle 

Subscripts 

c  Conical  flow 

e  Edge  of  boundary  layer  conditions 

w  wall  conditions 

s  stagnation  conditions 

o  pressure  induced 

t.c.  transverse -curvature  induced 
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Subscripts  (Cont'd) 

free  stream  conditions 
o  no  blowing 

Superscripts 

•  Reference  enthalpy  conditions 

-  Averaged  values  for  angle  of  attack  properties 
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INTRODUCTION 
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This  investigation  is  part  of  the  Nike-X  discrimination  studies 
which  are  aimed  at  developing  a  quantitative  evaluation  of  the 
observable  characteristics  of  a  re-entry  vehicle  and/or  decoy. 

The  drag  coefficient  of  a  vehicle  is  a  major  discriminant  because 
it  is  the  most  significant  parameter  affecting  the  vehicle's  velocity- 
altitude  signature. 

The  primary  purpose  of  this  investigation  is  to  develop  analytical 
expressions  for  evaluating  the  total  drag  coefficient  of  spheres 
and  cones.  The  accuracy  of  these  equations  is  to  be  determined  by 
comparing  the  results  with  wind-tunnel,  ballistic  range,  and  flight 
test  data.  Since  the  drag  equations  presented  in  this  report  are  to  be 
used  in  discrimination  studies,  it  is  necessary  that  the  independent 
variables  consist  only  of  free-stream  properties,  body  geometry,  and 
surface  material.  Consequently,  all  the  viscous  calculations  which  normally 
involve  .  flow  properties  at  the  edge  of  the  boundary  layer  have  to  be 
modified  in  order  to  transform  to  free-stream  conditions.  This  require¬ 
ment  generally  leads  to  the  employment  of  correlated  conical  shock 
relations,  perfect  gas  laws,  and  curve -fitting  techniques.  For  the 
inviscid  calculations  it  was  also  necessary  to  develop  curve-fits  which 
represent  the  results  of  rigorous  solutions  in  terms  of  the  desired 
parameters. 


In  the  employment  of  these  simplifying  techniques,  it  was  necessary 
to  maintain  an  accuracy  requirement  that  is  within  the  range  specified 
for  this  investigation;  i.  e. ,  10%  accuracy  from  sea  level  to  100,  000  ft, 
and  then  a  linear  variation  up  to  30%  accuracy  at  200,000  ft.  This 
accuracy  requirement  is  based  on  the  following  ranges  of  variables: 


Cone  half  angle 
Cone  length 
Cone  bluntness  ratio 
Sphere  diameter 
Wall  temperature 
Free-stream  Mach  number 
Altitude  or  approximate 
Free-stream  Reynolds  number 
Angle  of  attack 


8  to  15  degrees 
1  to  12  feet 
0  to  0.  3 

3.  75  to  12  inches 
1000  to  3000°  Rankine 
5  to  25 

0  to  200,  000  feet 
104  to  1010 

0  to  20  degrees  at  200,000  feet 
0  to  0  below  150,  000  feet 


The  accuracy  requirement  is  also  based  on  a  vehicle  that  is 

not  ablating  or  experiencing  mass  injection  into  the  boundary  layer. 

Since  most  re-entry  vehicles  will  have  some  type  o£  thermal 
protection  which  does  ablate,  the  accuracy  requirement  will  not  be 
met  if  mass  addition  effects  are  not  included.  Although  the  contractual 
requirements  for  this  investigation  did  not  include  mass  addition 
(blowing),  an  approximate  method  for  determining  the  reduction  in 
skin  friction  drag  due  to  blowing  is  provided.  This  approximate 
method  should  result  in  a  total  drag  coefficient  that  is  within  the  specified 
accuracy  range. 

Also  included  is  a  discussion  of  more  rigorous  methods  of  analysis 
which  may  be  employed  to  determine  the  effects  of  blowing  on  skin- 
friction  drag  and  induced  drag. 

For  the  case  of  pure  spheres,  the  total  drag  coefficient  is  obtained 
solely  from  an  empirical  correlation  of  experimental  data.  For  the 
case  of  cones,  each  individual  drag  contribution  is  separately  determined 
by  either  pure  theoretical  techniques  or  semi -empirical  methods.  The 
individual  components  are:  (1)  forebody  pressure  drag  (CDp)  -  the  drag 
due  to  direct  pressure  forces,  (2)  base  pressure  drag  (CDb)  *  tiie  drag 
due  to  a  pressure  in  the  base  region  which  is  less  than  ambient,  (3)  skin- 
friction  drag  (CDf)  -  the  drag  due  to  viscous  shear  forces,  laminar  and 
turbulent,  and  (4)  induced  drag  (Cj^)  -  the  increased  pressure  and  skin- 
friction  drag  due  to  a  thick  boundary  layer  (considered  for  laminar 
flow  only)  which  interacts  with  the  inviscid  flow. 

The  transition  from  laminar  to  turbulent  flow  has  been  investigated, 
and  predictions  are  made  for  determining  the  altitude  range  where 
this  phenomena  occurs.  F or  the  case  of  cones,  the  results  of  various 
flight  tests  have  been  employed  to  produce  a  transition  criterion. 
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II.  CONE  PRESSURE  DRAG 


A.  Sharp  Cone  Forebody  Pressure  Drag,  a  =  0 

The  forebody  pressure  drag  coefficient  is  defined  by  the  relation. 


2s 

ab 


/ * 

Jo 


For  a  sharp  cone  at  a-  0,  the  pressure  (pc )  is  constant  over  the  entire 
body,  and  the  forebody  drag  is  then  only  a  function  of  the  pressure  ratio 
( pc/p^  across  the  conical  shock.  Solutions  for  this  pressure  ratio  as  a 
function  of  Mach  number  and  cone  angle  have  been  tabulated  by  Kopal 
(reference  1}  for  a  ratio  of  specific  heats,  y  =  1.405.  These  results  have 
been  employed  directly  in  the  drag  equation  in  order  to  obtain  the  forebody 
drag  on  a  sharp  cone. 

The  results  have  also  been  compared  to  exact  conical  flow  solutions  em¬ 
ployed  at  Avco  RAD  for  the  case  of  a  real  gao  (Program  1427)  (reference 
2),  and  there  was  excellent  agreement. 

B.  Blunt  Cone  Forebody  Drag,  «  =  0 

For  a  blunt  cone  at  o=  0,  the  pressure  (pc)  in  the  forebody  drag  equation 
is  not  constant  with  respect  to  axis’  distance  (*).  However,  at  a  given 
value  of  x,  the  pressure  is  constant  in  the  transverse  direction  when  *  *  0. 
Therefore,  to  solve  the  forebody  drag  equation,  the  axial  pressure  dis¬ 
tribution  pc  *  fa[r(x)l  must  first  be  determined,  and  then  the  equation  can 
be  integrated  directly. 

The  pressure  distribution  on  a  blunt  body  is  most  accurately  determined  by 
the  method  of  characteristics,  which  is  explained  in  reference  3. and 
numerous  other  texts.  The  numerical  solution  of  this  method  may  be 
obtained  from  digital  computer  programs  developed  by  Avco  RAD.  To 
provide  the  initial  data  for  the  method  of  characteristics  (employed  in  the 
supersonic  region),  a  solution  for  the  transonic  flow  region  around  the 
blunt  nose  is  first  employed  (Program  1447  )  (reference  4).  The  initial 
data  is  then  used  in  the  axi -symmetric  flow  method  of  characteristics 
solution  (Program  596)  (reference  3). 

For  this  present  investigation  of  blunt  cones,  method  of  characteristics 
solutions  comparable  to  those  obtained  from  the  digital  computer  programs 
mentioned  above  were  employed. 
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These  solutions  were  originally  documented  by  Chuskin  and  Schulishnina 
in  the  form  of  tables;  their  results  were  translated,  edited,  and  re- 
uocumented  at  Avco  RAD  {reference  5  ).  These  solutions  were  obtained 
**v  the  method  of  characteristics,  employing  the  blunt  body  solution  of 
Beloisorkovskii  (reference  6)  for  the  initial  data. 

The  reasons  for  selecting  to  employ  these  solutions  in  this  investigation 
are:  1)  the  results  are  essentially  the  same  as  obtained  by  the  Avco  RAD 
methods,  2)  the  pressure  distributions  have  already  been  integrated  to 
yield  the  forebody  drag  coefficient,  and  3)  it  would  take  considerable 
computer  time  to  re-run  all  the  solutions  investigated  in  reference  6  . 

C.  Correlation  of  a  -  0  Solutions 

The  results  of  the  sharp  and  blunt  body  solutions  at  a  -  0  have  been  cor¬ 
related  as  a  function  of  0  .  and  RN/RB  .  (reference  7  ).  The  analysis 
consisted  of  curve-fitting  the  solutions  in  the  form. 


which  is  similar  to  the  form  of  the  Newtonian  drag  equation.  The  quantities 
C  ,  N,  K,  A,  and  B  represent  the  variations  from  Newtonian  theory  due  to 
Mach  number  effects.  Employing  polynomial  curve -fitting  techniques,  we 
obiain, 

C  -  1.944  +  1.872  MJ1  -  17.0C  M„~2  38.194  MJ"3 

N  »  1.931  +  0.8635  -  8.063  M~2  +  12.205  M„~3 

X  -  11.433  +  34.96  MJ*1  -  921.51  M„~2  26C7.3  *C3 

A  -  0.5359  +  0.09964  M„“l  +  10.769  MJ,"2  -  104.21  M„_3  +  209.43  MJ-4 

B  -  3.296  +  2.997  MJ*1  -  74.378  MJ*2  +  154.67  MJ"3 

The  results  of  these  curve -fits  are  compared  to  the  basic  method  of  char¬ 
acteristics  solutions  in  figures  1  through  4  for  Mach  numbers  of  5,  10,  17, 
and  25,  and  varying  bluntness  ratios.  The  agreement  is  good,  therefore 
satisfying  the  first  requirement  of  developing  an  analytical  expression  for 
forebody  drag  at  zero  angle  of  attack.  It  should  be  emphasized,  however, 
that  the  curve  fits  are  valid  only  within  the  specified  limits  of  M,  0 ,  and 

Rn^b  * 


Figure  2  DIGITAL  CURVE  -  FIT  COMPARISON  WITH  BANC  DATA;  FOWEBODy  DRAG 
VARIATION  WITH  CONE  ANGLE  AND  BLUNTNESS  PATIO  -  M  .  !0.' 


Figure  3  DIGITAL  CURVE  -  FIT  COMPARISON  WITH  BASIC  DATA;  FOREBODY  DRAG 
VARIATION  WITH  CONE  ANGLE  AND  BLUNTNESS  RATIO  -  M ,  17.0 
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In  order  to  account  for  angle  of  attack  effects  on  forebody  pressure  drag, 
the  following  methods  of  analyses  may  be  employed: 


1)  Sharp  Cone 

The  circumferential  pressure  distribution  around  the  yawed  body  (in 
a  plane  perpendicular  to  the  axis  of  symmetry)  may  be  obtained  from 
a  perturbation  solution  of  the  conical  flow  equations.  This  solution  is 
presently  programmed  for  machine  computation  (Program  1431) 
(reference  8  ).  The  resulting  pressure  distributions  are  then  inte¬ 
grated  over  the  body  to  yield  the  forebody  pressure  drag  coefficient. 


2)  Burnt  Cone 

The  circumferential  pressure  distribution  for  this  configuration  may 
be  obtained  from  a  recently  developed  method  of  integral  relations 
(reference  9  ).  Although  this  solution  has  been  programmed,  it  has 
not  been  completely  checkout  out;  therefore,  this  method  of  analysis 
has  not  been  employed  in  this  investigation. 


3)  Newtonian  Theory 


This  widely  accepted  theory  may  be  applied  to  both  sharp  and  blunt 
cones  at  angle  of  attack.  The  theory  may  be  used  directly  to  calculate 
forebody  pressure  drag,  or  it  may  be  used  to  generate  angle  of  attack 
scaling  factors  which  are  then  applied  to  the  a  =  0  solutions  presented 
in  the  previous  paragraphs. 


A  limited  comparison  has  been  made  between  methods  1  and  3,  and  the 
results  indicate  excellent  agreement  in  the  hypersonic  range  when  Newtonian 
theory  is  used  to  generate  scaling  factors;  i.e. , 


frp  )a  *  (CDp)a-0 


Newt 


# 


where  (Cq  \  is  obtained  f-c;n  the  correlations  presented  in  paragraph  C 
\  p/e  m  o 

of  this  section. 

The  scaling  factor  /Cr,  /Cn  \  is  obtained  from  Newtonian  theory 

\Ua  ua  m  0)  New* 

employing  digital  computer  program  814  (reference  101.  In  this  solution, 
the  Newtonian  pressure  distributions  at  angle  of  attack  are  integrated 
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circumferentially  and  longitudinally  along  an  axi -symmetric  body  of 

revolution.  The  program  computes  the  bounds  of  the  aerodynamic  shadow 

region,  and  calculates  aerodynamic  coefficients  for  angles  of  attack  from 

0  to  it  radians  at  any  input  interval  desired.  Some  results  of  this  program 

are  plotted  in  figures  5  through  7  as  the  drag  coefficient  CD  versus  the 

cone  angle  0  for  bluntness  ratios  R^/Rg  >  of  0,  0.  1,  0.2,  and  0.3  for 

angles  of  attack  of  0°,  10°,  and  20°.  Figure  8  presents  the  desired 

Newtonian  drag  ratio,  (Cn  /CD  )  ,  as  a  function  of  cone  angle, 

\  «  a  »  0'  Newt 

bluntness  ratio,  and  angle  of  attack. 


In  order  to  develop  an  analytical  expression  for  this  ratio,  the  results 
presented  in  figure  8  have  been  curve  fit,  employing  a  :  ecently  developed 
digital  computer  program.  Avco  RAD  program  1413  (reference  II  )  enables 
a  dependent  variable  to  be  fitted  as  a  function  of  polynomial  powers  of  up 
to  three  independent  variables.  The  program  uses  the  least-squares 
approach  to  data  fitting  where  the  independent  variable  are  all  presented 
as  polynomials  cf  a  specified  order. 


The  resulting  curve -fit  is  of  the  form: 

(C^CPo«a)wewt  “  ^Al  +  A2a  +  ^A4“  +  A5°  + 

+  '  f  Ag  (.  -*■  A<j  cr }  ! 


+  1  +  Ajj  a  A  Aj2  a2!  +  [Ajj  +  Aj^  a  v  Aj^  a  1 1  0 


+  lAjg  +  A17  a  *  Ajg  a2  ]  S2  1  — ■ 

KB 


+  UAjj  +  Ajq  a  +  Ajj  a2 1  +  (Aj?  A23  c  *  Ajg  fl2  1$ 


+  'A25  +  A26  0  +  A27  o2  1  ^  ^ 

where  the  numerical  values  of  the  coefficients  are  given  in  tat.e  I.  The 
accuracy  of  the- curve -fit  is  within  5  percent  of  the  actual  Newtonian 
values. 
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•  Burs  5  NEWTONIAN  PRESSURE  DRAG  COEFFICIENT,  0  ANGLE  OF  ATTACK 


NEWTONIAN  PRESSURE  DRAG  COEFFICIENT  i0!'  ANGi. 


TABLE  I 


COEFFICIENTS  FOR  (fcD  /CD  \  CURVE -FIT  SOLUTION 

V  a  a  -  d/Newt 


A1  * 

4.  54681695  E-01 

a2  » 

1.08259007  E-01 

a3  = 

7.  52646267  E-02 

a4  = 

9.92751992  E-02 

a5  = 

-  1.  33200192  E-02 

a6  “ 

-  9.  13259697  E-03 

A?  = 

-  4.  32304686  E-C3 

a8  = 

6. 08808160  E -04 

Ag  = 

2.95970723  E-04 

A10  = 

1.63346966  E  +  01 

All  = 

-  4.  50010377  £  +  00 

aI2  = 

-  1.07061625  E-01 

a13  = 

-  2.97740823  E  +  00 

A 14  = 

8.21443164  E-01 

A 1 5  = 

5.86190152  E-03 

A 16  = 

1.28462178  E-01 

A 17  = 

-  3.55578628  E-02 

A 18  = 

7. 0047 1497  E-05 

A 19  = 

-  5.  35168993  E  +  01 

A20  * 

1.51758946  E  +  01 

A2I  = 

-  4.  10820156  E-01 

a22  = 

9.74715936  E  +  00 

a23  ' 

-  2.77304727  E  +  00 

a24  = 

8.  15759611  E-02 

A25  =  -  4.20040721  E -01 
A26  =  1.09716924E-01 

A27  =  -  3.72121710  E-03 


in  order  to  confirm  the  validity  of  Newtonian  scaling  laws  at  angle  of  attack 
a  typical  comparison  with  experimental  data  is  presented  in  figures  9  and  10 
The  comparison  consists  of  plots  of  the  pressure  coefficient  (Cp)  versus 
meridian  angle  (d)  for  sharp  10s  and  blunt  15®  cones  at  zero  and  10°  angle 
of  attack.  The  Newtonian  pressure  coefficient,  Cp  was  obtained  from  the 
relation, 


P  -  P» 


for  the  sharp  cone  case,  the  result  of  integrating  these  pressure  coefficients 

yielded  a  ratio  /Cn  ^  that  was  only  2  percent  higher  than  the 
\  ua  ua  -  0/Newt 

correspondhg  ratio  for  the  experimental  data.  For  the  blunt  cone  case 
(figure  10),  the  Newtonian  drag  ratio  is  3  percent  higher  than  the 

experimental  drag  ratio. 


The  data  of  figures  9  and  10  in  conjunction  with  other  similar  comparisons 
(reference  12,  13)  have  been  correlated  to  determine  the  error  incurred  in 
employing  Newtonian  theory  for  pressure  distributions.  The  results  of  this 
correlation  indicate  that  at  relatively  low  Mach  numbers  and  small  angles 
of  attack,  Newtonian  flow  under -estimates  the  pressure  distribution  on  the 
windward  side  of  the  body.  As  both  the  Mach  number  and  angle  of  attack 
are  increased,  the  comparisons  show  improvement.  The  Newtonian 
assumption  of  a, zero  pressure  coefficient  for  the  leeward  side  in  the  shadow’ 
region  tends  to  result  in  overestimation  of  the  leeward  pressures  by  as 
much  as  a  factor  of  two  in  some  cases.  However,  the  leeward  pressures 
at  angle  of  attack  are  of  course  very  much  smaller  than  the  windward 
pressures,  and  accordingly,  the  effect  on  vehicle  drag  of  relatively  large 
errors  in  leeward  pressure  will  not  be  substantial. 


In  order  to  present  some  quantitative  comparisons,  the  percentage  error 

in  C  defined  bv 
P 


has  been  computed  for  the  windward  meridian.  Figure  i  1  shows  a  plot  of 
this  percentage  error  as  a  function  of  the  similarity  parameter  K  defined 
as 


K  «  2  tan  (6  +  a) 


(7) 


-16- 


371# 


Also  included  in  figure  1 1  are  comparisons  of  Newtonian  theory  with 
analytical  values  computed  from  exact  conical  flow  solutions,  at  «  -  0 
as  given  in  refe  *ence  14.  The  percentage  error  for  this  latter  information 
is  defined  as: 


Figure  11  shows  that  Newtonian  theory  underpredicts  windward  pressures 
at  the  lower  K  values.  The  underprediction  is  about  10  percent  for  a  K 
value  of  about  3.  0.  For  a  10  degree  cone,  at  a  5  degree  angle  of  attack, 
for  example,  a  K  value  of  3.0  would  occur  at  a  Mach  number  of  5.6.  At 
K  values  of  the  order  of  10.  0,  the  Newtonian  theory  overpredicts  the 
pressure  coefficient  by  about  3  percent.  For  the  same  10  degree  cone  at 
5  degree  angle  of  attack,  a  K  value  of  10. 0  would  be  attained  at.  a  Mach 
number  of  18.7.  Thus  it  can  be  deduced  that  for  moderate  cone  angles  and 
even  relatively  small  angles  of  attack,  the  Newtonian  theory  tends  to  be  a 
reasonable  approximation  over  the  hypersonic  flow  regime. 


For  an  8  degree  cone  at  zero  angle  of  attack  at  Mach  5.  0,  which  represents 
the  lower  limit  of  the  K  parameter  in  the  current  study,  and  therefore  the 
maximum  error  in  Newtonian  flow,  the  percentage  error  is  seen  to  be  about 
24  percent.  However,  this  is  the  error  incurred  when  Newtonian  theory 
is  employed  directly  in  the  absolute  sense.  In  this  investigation,  Newtonian 
theory  is  only  used  to  scale  the  zero  angle  of  attack  solutions  which  are 
obtained  by  more  accurate  methods.  Consequently,  for  an  8  degree  cone  at 
»!.,=  5,  the  error  will  be  much  less  than  indicated  by  figure  11,  It  is  also 
shown  in  this  graph  that  as  any  of  the  parameters,  Mach  number,  cone  an 
angle,  or  angle  of  attack,  one  increased,  the  percentage  error  incurred 
from  Newtonian  theory  decreases. 

E.  Base  Pressure  Drag 

The  basic  data  used  to  calculate  base  drag  coefficients  were  obtained  from 
flight  test  pressure  measurements  reported  in  reference  15.  This  flight 
test  data  has  been  verified  by  ground  test  data  reported  in  references  16 
and  17.  Investigation  of  the  data  shows  that  the  majcr  parameters  in¬ 
fluencing  base  pressure  drag  are' Mach  number  and  cone  angle.  The  effects 
of  Reynolds  number,  wall  temperature,  bluntness  and  ablation  rates  did  not 
significantly  effect  the  average  pressure  level  in  the  separated  flow  region 
of  the  base. 


ACCURACY  OF  PROCEDURE  FOR  CALCULATING  SURFACE-PRESSURE  COEF 
»NE  CORRELATED  BY  MEANS  OF  HYPERSONIC  SIMILARITY  PARAMETER 


To  correlate  the  effects  of  Mach  number,  the  flight  test  results  obtained 
from  reference  15  have  been  used  to  define  a  reference  base  drag  co¬ 
efficient  (Cp^)^  •  *or  a  cone  angle,  0  =  0  degree.  A  curve -fit  of  the 

variation  of  (Cn  )  .  with  Mach  number  has  been  reported  in  reference  7 
.  .  UB  tel 

to  be 

/c  \  „  e(-1.706  -  0.33^  +  0.00566  mJ) 

\  DB/  ref 

This  variation  is  presented  in  figure  17*  together  with  the  limiting  vacuum 
condition. 

To  correlate  the  effects  of  cone  angle,  the  results  of  reference  12  show 
that  Cr.  •»  (1  +  sin  d ) . 

The  final  equation  for  the  base  drag  coefficient  i3  then, 


% 


(1+*“MC°B)cef 


„  (-1.706  -  0.33  M  +  0.00566  ) 

(1  +  sin  0)  e 


* 


s 


MACH  NUM&SHH  ,  M  ^ 


igvre  12  REFERENCE  BASE  -  DRAG  COEFFICIENT  (CDB)R£F  F 


m.  CONS  SKIN-FAICTION  DRAG 


The  skin-friction  drag  coefficient  is  defined  by  the  relation. 


'i( 


Cf  cos  6  rdx 


where  C{  is  the  local  skin-friction  coefficient  defined  by  the  relation. 


2 

—  p_  u_ 

2  "  “ 

The  following  two  sections  are  devoted  to  the  derivation  of  the  equations  for 

CD  in  both  laminar  and  turbulent  flow  conditions  for  a  sharp  cone  at  a  a  0, 

f 

A.  Sharp  Cone,  Laminar  Flow,  a  s  0 

To  compute  the  local  skin-friction  coefficient,  we  employ  the  Blasius  flat 
plate  incompressible  solution  modified  for:  1)  conical  flow  by  the  Mangier 
transformation  (reference  18):  2)  compressibility  by  Eckerts  reference 
enthalpy  method  (reference  19);  and  3)  conversion  from  conditions  at  the 
edge  of  the  boundary  layer  to  free-stream  conditions;  i. e. , 

P e  ue2 

Cf-  • 

The  modified  Blasius  equation  is  then. 


P *  m"  (p*  ***  \ 
P*  P*  \Pm  uj/ 


0.5  „  1.5  „  0.5 


.0.5  .0.5 


.0.5  .0.5 


Setting  =— 


Pe  V  Pt  P  P 
Pm  Re0'3  Pe°"3 

and  re-arranging,  we  obtain. 


Pm  «*• 
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Employing  the  curve -fit  results  of  reference  20 


Since  p*  *  pe  ,  and- 


l  ,*  /T.\ 

17 *X  ^  (*••«=>!»< S.  *SJ 

•,  •  Tcfer  fer  ©■“" 


©■ 


where 


h_ 

V 


0.5  4.  0.5 


+  0.22  r 


(51)  < 


\-i 


",io‘'  (£)  ’  (r)  •  “d  (r) ; 

_ !_•  l ' .  tL a. 3  ^  _  / 


*■  0.0374  Mez  msnming  - 


1.4,  i  m  0.85 


s.  ■  s.  ■  s. 


are  obtained  from  conical  flow  re  stilts 


which  have  been  correlated  as  a  function  of  the  hypersonic  similarity  param¬ 
eter,  Mm  sin0.  These  correlations  are  valid  for  the  range  of  this  investiga¬ 


tion  which  is  0.  7  —  tin 6  -  6.  5.  This  range  is  specified  by  an  8*  cone 


at  M„  5  5  as  the  lower  limit,  and  a  15*  cone  at  Mm=  25  as  the  upper  limit. 

The  velocity  ratio  (ue  /uw)  is  obtained  from  the  conical  flow  results  presented 
in  reference  21,  A  curve-fit  of  these  results  yield: 


i.4 


sin  0) K9 


0.5 


The  pressure  ratio  (pe  lpm )  is  obtained  from  the  curve-fit  relation  presented 
in  reference  22,  ~  X 


re 


P- 


-1+2.8  M„2  sin*  0 


*2.5  +  1  sta  0* 

1  +  16  sin0 


The  temperature  ratio  (Te/T„)  is  obtained  from  the  conical  flow  solutions 
presented  in  reference  1.  A  curve-fit  of  these  results  yield: 


-  -  1  ♦  0.0966  M  sin  6  +  0.2267  (M  sin  9) 

T« 

The  local  Mach  number  is  then  computed  from, 

-0.5 


"* '  “•  fe)  © 


Since  the  local  skin-friction  coefficient  is  now  defined  in  terms  of  free- 
stream  conditions,  wall  temperature,  and  cone  angle,  we  may  now  integrate 
over  the  cone  surface  to  obtain  the  skin-friction  drag.  Therefore, 


B'^il 


Cf  cos  9  (2 or)  dx 


Since 


and  t  =  x  sin  9 


c«.  ■  >•  (r) 

Cp^  .  — ~ ( Cj  1,  _  '  co.  6  .in  8  f  .1. 


4  cos  9  sia  9 

•7  “77“  L 


Since 


Rb  -  L2  sin2  9 


and 


we  cbtain  the  final  result  for  laminar  flow. 


1.53  U  ^  °'5  ^  ^  !85  ^  -O'58’ 
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cot  9 


The  solution  of  this  equation  is  presented  in  figure  l*  for  8*  and  15*  cones. 
B.  Sharp  Cone,  Turbulent  Fow,  a  *  0 


We  shall  employ  the  Blasius  flat  plate  incompressible  solution  moified  by 
the  conical  flow  transformation  factor  reported  in  reference  23,  modified 
for  compressibility  by  the  reference  enthalpy  method  (reference  19)  and 
converted  to  free -stream  conditions.  The  modified  equation  is  then. 


and  re-arranging. 


tU3S„3  . 

X 


0.2 


Again  Employing  the  cui  t>e -fit  results  of  reference  20,  modified  slightly  in 
order  fij  adequately  cb^er  the  range  of  h  ,  he  ,  and  h„  encountered  in  this 
study,  we }. 


^y-2  ^  0  58  j°a 


Assuming  Cp 


=  Cp  a  Cp  ,  and  setting  y  *  1,  4  and  r  -  0.  38 
*  003“"' 


The  local  skin-friction  equation  in  terms  of  free -stream  conditions,  wall 
temperature,  and  cone  angle  is  now. 


-“ferw’er 


X 

where  the  "e"  conditions  are  determined  from  the  same  equations  specified 
in  the  previous  paragraph. 


To  obtain  the  skin-friction  drag  coefficient,  the  local  coefficients  are  inte¬ 
grated  over  the  cone  surface  with  the  condition  that 

c'-  -  *0-1  w  • 
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FRICTION  DRAG  COEFFICIENT,  SHARP  CONE,  LAMINAR  FLOW 


Therefore, 


CDf  “  ~  (CfJx-L  L°-2  cos  6  si: 


/ 


0.8  i 

x  dx 


Since 


10  cos  0  sin  6  - 

- (Cf  )  »  -  L2 

RB 


Rn  ■  L2  sin2  6 


and 


(Cf  ). 


'f„x  *  L 


0.0698  /Ue\18  /Pe\°‘8  /Te 

(R  >°-2  yj  \pj  \T- 


^  -0.58  /  P-58 

I  \s) 


we  obtain  the  final  result  for  turbulent  flow, 


Tie  solution  of  this  equation  is  presented  in  figure  14.  for  8°  and  15°  cones. 
C.  Bluntne as  Effects 


The  effects  of  nose  bluntneos  tend  to  reduce  the  skin-friction  drag  compared 
to  the  sharp  body  case.  This  reduction  is  due  primarily  to  the  loss  in 
dynamic  pressure  encountered  by  the  free  stream  as  it  passes  through  a 
strong  curved  shock.  An  additional  reduction  in  skin-friction  is  caused  by 
the  decreased  cone  pressure  resulting  from  overexpansion  of  the  flow. 

To  completely  evaluate  the  effects  of  nose  bluntness,  the  following  methods 
of  analysis  have  been  employed: 

a.  Method  of  characteristics  (reference  3  )  (Avco  Computer  Program 
596).  This  method  was  used  to  obtain  the  pressure  distribution  and 
shock  shape. 

b.  Blunt  Body  Viscid-Inviscid  Solution  (references  24  and  25)  (Avco 
Computer  Program  1115B).  This  method  was  used  to  compute  the 
distribution  of  flow  properties  along  the  edge  of  the  boundary  layer  and 
the  resulting  local  skin-friction  coefficients. 
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In  the  employment  o£  the  latter  method  o£  analysis,  the  known  parameters 
are  free  stream  conditions,  body  geometry,  surface  temperature,  surface 
pressure,  and  shock  shape.  The  additional  fluid  property  that  can  be  eval¬ 
uated  along  the  edge  of  the  boundary  layer  is  entropy.  This  is  determined 
from  an  iterative  solution  in  which  the  free -stream  mass  flow  is  equated 
to  the  mass  flow  in  the  boundary  layer.  The  analytical  model  used  to  per- 
form  tills  computa  tion  is  shown  in  figure  15. 


To  obtain  the  flow  properties  at  the  edge  of  the  boundary  layer  at  Station 
"X",  a  value  for  the  entropy  is  initially  assumed.  Since  the  pressure  is 
known,  all  other  properties  may  be  computed.  The  mass  flow  in  the  bound¬ 
ary  layer  is  then  computed  by  employing  the  correlations  reported  in  refer¬ 
ence  24  which  permit  an  approximate  calculation  of  the  mass  flow  as  a  func¬ 
tion  of  boundary  layer  edge  conditions  and  surface  temperature.  Employing 
the  continuity  equation  between  the  boundary  layer  at  Station"X"  and  the  free 
stream,  the  bounding  stream-tube  radius,  R„,  is  determined.  From  the 
known  shock  angle  corresponding  to  Rm,  the  entropy  change  across  the  shock 
can  be  computed.  The  resulting  entropy  is  compared  to  the  assumed  value 
and  iteration  is  made  until  convergence  is  achieved.  The  final  fluid  prop¬ 
erties  at  the  edge  of  the  boundary  layer  are  then  used  in  the  modified' 
Blasius  equation  in  order  to  compute  the  local  skin-friction.  The  same 
procedure  io  then  employed  at  the  next  station  along  the  body.  After  all 
stations  are  computed,  the  local  skin-friction  coefficients  are  integrated 
over  the  surface  of  the  body  in  order  to  obtain  the  skin-friction  drag  coef¬ 
ficient. 


The  foregoing  methods  of  analysis  have  been  employed  in  a  parametric 
study  in  order  to  develop  simplified  correlations.  The  major  problem  en¬ 
countered  in  this  study  is  the  number  of  variables  affecting  blunt  body 
viscous  drag;  they  are;  Rjt/P.g  ,  R«L,  M«. ,  0,  and  Therefore  a  prelim¬ 
inary  investigation  was  made  of  the  relative  significance  of  each  parameter 
in  order  to  determine  which  parameters  were  the  most  critical  and  which 
ones  could  be  neglected. 


For  laminar  flow  it  was  determined  that  the  critical  parameters  were  rj/rb> 
and  R„. .  Therefore,  the  study  was  conducted  for  fixed  values  of  Kb.  s  22, 

0  =  10  ,  and  Tw  =  3,G00*R.  The  results  are  presented  in  figure  16  where 
the  blunt  body  skin-friction  drag  is  normalized  by  the  corresponding  sharp 
cone  value.  The  characteristic  length  in  the  Reynolds  number  is  the  slant 
length  (L)  of  the  sharp  cone.  When  the  known  quantities  are  cone  angle, 
bluntness  ratio,  and  blunt  body  axial  length,  the  slant  length  is  determined 
from  the  relation. 


la 

L  «=  -  ^  • 

rn 

cos  S  —  -  (1  -  tin  6) 

RB 
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Figuro  15  ANALYTICAL  MODEL  FOR  BLUNT  BODY 
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FRICTION  DRAG  RATIO,  LAMINAR  FLOW 


The  values  of  R„L  shown  in  figure  16  represent  many  combinations  of 

altitude  and  body  length,  and  in  all  cases  the  correlation  was  quite  satis¬ 
factory.  To  further  confirm  this  correlation,  a  few  cases  were  investigated 
with  the  velocity  reduced  by  one -half.  The  resulting  correlation  was  very 
good  up  to  =  10^*;  at  higher  Reynolds  numbers  there  wao  some  deviation 

due  to  Mach  number  effects. 

The  results  presented  in  figure  16  show  some  interesting  trends  which 
are  summarized  as  follows: 

(1)  At  low  values  of  R„j  (<  10^),  the  bluntness  effects  are  negligible 
due  to  the  thick  boundary  layer. 

(2)  At  high  values  of  (>  10^),  the  boundary  layer  is  so  thin  that  the 
streamline  intersecting  the  edge  of  the  boundary  layer  at  *  =  L  has 
passed  through  a  region  of  the  shock  that  is  essentially  normal.  There¬ 
fore,  the  minimum  value  of  the  drag  ratio  corresponds  to  the  case 
where  the  entropy  along  the  boundary  streamline  equals  the  entropy 
behind  a  normal  shock. 

For  turbulent  flow,  the  parametric  study  indicated  that  the  effect  of  R~L 
was  not  critical  in  the  high  Reynolds  number  range  {>  105).  The  dominating 
effects  were  found  to  be  the  bluntness  ratio  and  free-stream  Mach  number. 
The  other  effects  of  cone  angle  and  wail  temperature  were  assumed  to  be 
negligible  in  the  case  of  laminar  flow.  The  resulting  variation  of  the 
skin-friction  drag  ratio  is  presented  in  figure  17.  These  results  show 
that  the  ratio  steadily  decreases  up  to  a  bluntness  ratio  of  about  0.  3,  and 
then  gradually  approaches  its  asymptotic  limit,  which  is  the  case  of  normal 
shock  entropy  at  the  edge  of  the  boundary  at  *  =  L, 

The  results  presented  in  figures  16  and  17  have  been  curve-fit  to  yield 
the  following  equations: 

Laminar  Flow 
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with  the  conditions; 


»f  A  >  1.0,  set  equal  to  1.  0 

if  R  <  R  •  ;  set  equal  to  R  . 

rain  *  cun 


where  R  .  =1.0-  6.5076 
mm 


A  /  rn\2  /*nV 

-1+  30.46  -  »  46.  85 - 

5/  \Hj  \RB  / 


Turbulent  Flow 


(C°f)  Blunt 


fDf)si 


1.0-  I  (0.8  +  0.052  #4  ) 


harp 


D.  Angle  of  Attack  Effects 


There  is  presently  no  reliaqle  solution  available  for  the  calculation  of  skin- 
friction  drag  on  yawed  axi-s/mmetric  bodies.  According  to  the  results  of 
reference  26,  the  cross -flow  effects  on  the  wi.jdward  and  side  meridians 
are  relatively  small  up  to  yaw  angles  equal  to  the  cone  angle.  This  implies 
that  for  a  -  0,  the  Blsaius  solution  can  be  employed  with  tangent-cone  ap¬ 
proximations  for  evaluating  the  fluid  properties.  Since  a  requirement  of 
this  study  includes  the  evaluation  of  skin -friction  drag  on  sharp  and  blunt 
cones  up  to  20*  angle  of  attack,  the  assumption  will  have  to  be  made  that 
the  modified  Blasius  solution  is  still  valid  over  this  entire  range  when  the 
fluid  properties  are  evaluated  using  tangent*>cone  approximations  for  the 
windward,  leeward,  arid  side  meridians.  Therefore,  for  the  case  of  a 


sharp  cone  in  laminar  flow, 


1  +  0.7  uj 


j  sin^  ( 0  +  a) 

I 


2.5  -  8  sin  (d?  +  a) 
1  +  16  Mw  sin  {{?  ♦  u) 


sin^  (6  —  a) 


2.5  tin  (0  -  a) 

1  +  16  sin  (<?  —  a) 


+  2  sin^  6 


2.5  +  8  M„  sin  0 
1  +  \6Umiiad 
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i  u 


f 


— •  l  +  0.02413  M  i*in  (Q  *  a)  +  sic  (0  —  a)  ♦  2  *in  0] 

T 

0* 

+  0.0567  uj  (*in^  (0  *■  a)  +  *in^  (6  —  a)  +  2  6l 

T  ‘  o  wo-5  ©fe)  *  0  0174  5'2 

wi'.h  the  condition  that, 

if  (0  -a)  <  0,  set  equal  tc  zei  ... 

For  the  case  of  turbulent  flow,  this  same  procedure  for  evaluating  the  aver¬ 
age  fluid  properties  is  utilized  in  the  zero  angle  of  attack  equation.  There¬ 
fore, 


[f°<>  J 


0.0776 
r~l  0,2 


(r  terr'er 


For  the  case  of  the  blunt  cone,  it  is  assumed  that  the  skin  friction  drag 
ratios  presented  in  figures  16  ar.d  17  are  still  applicable. 

E.  Mass  Addition  Effects 

The  effect  of  mass  addition  due  to  ablation  tends  to  reduce  the  skin-friction 
drag.  To  rigorously  account  for  this  effect,  the  following  analyses  should 
be  employed: 

1)  Evaluation  of  thermodynamic  and  transport  properties  of  the  in¬ 
jected  gas,  employing  Avco  RAD  digital  computer  programs  1291  and 
1619B.  (References  27  and  28). 

2)  Using  the  thermodynamic  and  transport  properties  in  the  solution 
of  the  boundary  layer  equations  for  both  laminar  and  turbulent  flow 
(Avco  RAD  digital  computer  programs  1475  and  1356,  respectively, 
references  29  and  30). 

i)  Correlation  of  the  skin-friction  reduction  ratio  (Cf/Cf  )  as  a  func- 

2  a  ® 

tion  of  the  blowing  parameter, - - 

ue  cf_ 
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Since  the  original  requirements  for  this  drag  study  did  not  include  mass 
addition  effects,  the  .rigorous  methods  of  analysis  have  not  been  employed. 
However,  provided  herein  is  a  simplified  approximate  method  which  yields 
reasonable  results  for  most  heat  shield  materials.  This  method,  which 
has  been  reported  in  reference  31,  employs  the  following  assumptions: 

a)  The  mass  injection  rate  (a)  can  be  computed  from  an  effective  heat 
of  ablation  (q*  }  and  the  heat  flux  to  a  non-ablating  wall  ).  Therefore, 

. 

m  «  — — 

q 

b)  The  heat  of  ablation  can  be  obtained  frem  the  relation, 

q*  -  Cp  AT  +  H,  +  <j(h,  -  hw) 

(total  energy  absorbed)  (energy  lost  by  transpiration) 

«  Hr  +  >?(hs  -  N,) 

c)  Reynolds  analogy  is  valid;  i.e.,  C(  f  C(  -  q/^ 

O 

The  energy  balance  required  at  the  wall  is  then, 
q  «  m  Rj. 

%  %  «T 

•  — —  H-*-  *  . - 

^  Ht  +  rj  (hs  -  b*) 

Solving  for  q/qp  ,  and  setting  equal  to  Cf/C{  ,  we  obtain, 

Cf  %  _ 

ci0  “t  +  i  <N  -  N,) 


This  skin-friction  reduction  ratio  due  to  mass  addition  effects  is  presented 
in  figure  18  for  laminar  and  turbulent  flow  for  varinus  heat  shield  mate¬ 
rials.  The  ablation  properties  used  in  this  approximate  method  are: 


Material 

t.CR) 

Hw/RT, 

Hy(Bm/lb) 

^Turb 

7L*o 

LT0 

1370 

9.9 

745 

0.41 

0.57 

Teflon 

1800 

13.4 

800 

0.255 

0.43 

OTWR 

4850 

42.0 

3200 

0.32 

0.50 
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IV  CONE  INDUCED  DRAG 


The  induced  drag  resulting  from  a  thick  boundary  layer  around  a  slender  body 
has  been  considered  by  Probstein  (reference  32)  and  others  for  the  case  of 
laminar  flow  with  no  mass  injection.  It  has  been  shown  by  Probstein  that  for 
sharp  slender  cones,  the  induced  effects  may  be  computed  from  "weak  inter¬ 
action"  considerations;  i. e.,  the  induced  pressure  gradient  does  not  significantly 
effect  the  boundary  layer  growth.  For  slightly  blunted  cones,  there  is  a  region 
of  "strong  interaction"  starting  immediately  att  of  the  stagnation  point  and  ex¬ 
tending  a  few  nose  radii  downstream.  However,  for  purposes  of  computing  in¬ 
duced  drag,  this  region  has  a  relatively  small  effect,  ami  weak  interaction  as¬ 
sumptions  may  be  employed  throughout. 

The  total  induced  drag  can  be  considered  to  be  the  summation  of  three  effects 
as  follows: 

1.  Induced  Pressure  Drag 

This  effect  is  due  to  a  new  apparent  body  shape  represented  by  the  displace¬ 
ment  thickness  profile.  Employing  the  tangent-cone  assumption,  the  new 

dS  » 

effective  cone  angle,  =  &  + - .  Since  d 5 /dz  decreases  as  *  increases, 

d* 

the  induced  pressure  is  greatest  near  the  leading  edge,  and  decays  at  <* 
rate  which  is  proportional  to  yfV. 

2.  Pressure  Induced  Skin  Friction  Drag 

This  effect  is  due  to  the  increase  in  local  skin  friction  as  a  result  of  in¬ 
creased  local  pressure.  As  in  the  case  of  induced  pressure,  this  effect  is 
greatest  near  the  leading  edge. 

3.  Transverse -Curvature  Induced  Skin  Friction  Drag 

This  effect  arises  when  the  boundary  layer  thickness  is  not  negligible  com¬ 
pared  to  the  body  radius  of  curvature.  The  negligibly  thin  boundary  layer 
assumption  was  made  by  Mangier  (reference  18)when  he  transformed  the 
boundary  layer  equations  for  an  axi -symmetric  body  to  a  flat  plate,  there¬ 
fore,  the  transverse -curvature  correction  is  used  to  modify  the  well-known 
Mangier  transformation  factor. 

The  following  paragraphs  present  the  derivations  employed  to  compute  in¬ 
duced  drag  in  terms  of  free  stream  conditions.  Results  are  presented  for 
laminar  flow  only,  since  induced  drag  is  negligible  in  turbulent  flow  when 
6  -  8*. 


A.  Induced  Pressure  Drag,  Sha  -p  Cone,  «  =  0 

Employing  Probstein's  equation  for  the  local  increased  pressure  on  a  ahar 
cone  (neglecting  second  order  terms  of  the  series  expansion), 

AP 

~  -  Fj(K)d.  Xe 
where: 


VCA 
V  V 


3 


with 


Ce  * 


0.968 

U" 


T_ 


♦  0.058  (For  y  «  1.4  And  Pr  «  0.72) 


Fj  OC) 


1  <Kpe/p,J 

/M 

_M~  i0 

»-9e  VV 

which  is  a  function  of  M„sin0  and  y.  A  curve-fit  of  Fj(K)  for  y  =  1.  4  yields, 
Fj  (JC)  »  0.9  -  0.119  -  0.0108  sin  8)2 

In  order  to  obtain  the  viscous  interaction  parameter  xcin  terms  of  free- 

JEL 


st  ' ,m  conditions,  \m  - 


mA  we  write: 


X- 

whe  re : 


X\0'6 


(5-)  '  & 

£  -  ffl 

i.iJ*  .pi) 

C-  *e  T»  V-/ 

Me  ue  /Te\-« 

VTJ 


e'0J 


(This  viscosity -tetapersture  power  ls» 
is  assumed  throughout ) 
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Substituting, 


The  ratios,  pe  /p^  ,  Te  /Tm,  and  uf  /uw  are  a  function  of  M„.  sin  6  and  are  oh 
tained  from  the  relations  given  in  section  in. 


Therefore,  the  equation  for  local  induced  pressure  with  all  quantities  in 
terms  of  free -stream  conditions  is: 


To  obtain  the  induced  pressure  drag,  it  is  then  necessary  to  integrate  this 
equation  over  the  cone  surface.  To  account  for  nose  bluntness  effects, 
some  approximations  must  be  made  in  order  to  integrate  the  equation  in 
closed  form.  (1)  The  value  of  the  induced  pressure  at  the  sphere-cone 
.angent  point  is  assumed  to  equal  the  vahis  for  the  corresponding  sharp 
cone  case  for  the  local  Mach  number  and  temperature  associated  with  the 
blunt  body,  (2)  The  local  flow  properties  are  assumed  to  be  invariant  with 
i,  and,  (3)  the  induced  pressure  on  the  spherical  nose  is  a  negligible  amount. 

The  body  geometry  employed,  for  th<  induced  drag  calculations  is  shown 
in  figure  19  • 

The  integration  will  now  be  performed  in  such  a  manner  that  when  Rj^Rg  =  0, 
the  solution  degenerates  to  the  sharp  body  case. 
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R^  cos  6 
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(K  x“0,5)  pe  t  df 


R^  cos  6 
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Pe 
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3  UP  RB2 


4  /  sin  8\0'5 

-T%  *Xt) 


-  VJ, 

7  J 


1  -  £ -  CO*  C7, 

RB 


Since - r=L»  w*  obtain 

sin  0 


ACr 


Crj  K 


/  rN  \  5/2 

-  W  "'1 


-0.5 


B.  fod»c«d  Skin-Friction  Drag.  Sharp  Cone.  «  =_j> 

Employing  Prob.teln'.  fir.,  order  elution  for  the  loci  Indued  .kin- 
_ -rn-ffi rinnt.  and  a*suminR  y-  1.4,  we  obtain, 

Fl(K) 


friction  coefficient,  and  itsuBiing  y 

r  VT 


<Acf>t> 


-f 


-  0.82}  +  0.52« 


VT- 


♦  0.458  M  ' 


Me' 


-43- 


where  d„,  Fl{K),  and  Xt  are  evaluated  in  the  acme  manner  aa  for  the  caae 
of  induced  pressure. 


Therefore, 


Ff) 


p  . 


where: 


[-0.823.0.524  — 

4-Fl00  uj  /T^  r«\"°'3/MW 

Vv  \pj  vj  W 


+  0.430  Me^ 


Employing  the  same  blunt  body  approximations  as  previously  discussed, 
we  may  now  integrate  in  closed  form: 
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Since 


//f 
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we  obtain: 
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z  «  i 


-44- 


for  the  non-interacting  *harp  cone  case,  it  may  be  shown  that  (aee  aeetion 

111): 


*D{  “  T  (C(>  * 


I  0.5 


CM  8 
*%a  6 


-0.5 


or 


(Cf>  .  -  — 

*  x  -  i  4 


3  3io«  ,  r  0.5  [0.5 
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CD, 


Substituting  in  the  induced  drag  equation, 

,  3  v'  r  ,ip2  8  (L1-5  -  i  L0-5 ) 

(ACd)  .  „  k  cL(  •  (L 


B 


Since 


N 


R  2  «  L2  «in2  6  and  l  -  L  co*  0  , 


we  obtain 


ho, 


c.  Traverse  Curvature  Induced  Skin-Friction  Drag,  Sharp  Cone,  a  -  0 

Employing  Probstein's  first  order  solution  for  the  local  induced  transverse 
curvature  effect,  and  again  assuming  y  s  1.4,  we  obtain 

Xe 


<ACf>t.c. 


21 

0.517  +  0.913  -r-=—  +  0.0484  ' 
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yT  tan  8 


-  K "  x 
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whe  re : 
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The  integration  is  performed  in  exactly  the  same  manner  as  the  pressure 
induced  skin-friction  drag.  Therefore, 


K)t.c.  -  4  L_0-5 


The  results  of  some  induced  drag  calculations  are  presented  in  figures 
20  to  23  for  sharp  8C  and  15°  cones. 

These  results  consist  of  induced  pressure  drag  (ACD  )  and  the  total  induced 
skin-friction  drag,  (ACD{ )  +  (ACDf)  .  p 

The  calculations  indicate  the  following: 

ACDf  ACd 

a.  The  ratio— —  is  significantly  greater  than  the  ratio— — . 

S  «  Cdp 

b.  The  more  slender  cone  experiences  much  higher  induced  drag  than 
the  steep  angle  cone. 

c.  Induced  drag  increases  with  wall  temperature. 

d.  The  effect  of  increasing  free-stream  Mach  number  generally  in¬ 
creases  the  induced  drag  (however,  there  is  a  reversal  effect  at  high 
values  of  wall  temperature). 

The  total  induced  drag  coefficient  is  then, 


C°i  '  +  (ACof)P  +  (ACDf)  t.c. 


D.  Bluntness  Effects 

For  the  case  of  the  blunt  cone,  the  flow  conditions  at  the  edge  of  the  bound¬ 
ary  layer  vary  with  x  for  two  reasons;  (I)  the  pressure  distribution  is  not 
s  pnstant  as  in  the  sharp  body  case,  and  (2)  the  curved  shock  produces  an 
er  ‘.ropy  gradient  along  the  edge  of  the  boundary  layer. 

Since  the  integration  scheme  presented  in  the  previous  paragraphs  is  valid 
only  for  constant  flow  properties,  several  assumptions  must  be  made  in 
evaluating  effective  constant  conditions.  This  approach  is  very  approximate, 
but  it  is  justified  on  the  grounds  that  for  a  blunt  body  the  induced  drag  is  a 
very  small  portion  of  the  total  drag;  in  fact,  it  is  most  often  neglected  with¬ 
out  any  consequences. 
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In  this  approach  the  major  assumptions  employed  are: 

a.  The  surface  pressure  ratio  for  the  non-interaction  case  (pe  /p^ ) 
is  the4  same  as  for  the  corresponding  sharp  cone  between  the  limits, 
f  <  «  <L. 

b.  The  assumed  constant  values  of  temperature  and  local  Mach  num¬ 
ber  are  to  be  evaluated  at  the  end  of  the  cone  (  x  =  L).  This  assumption 
was  made  in  reference  33  and  was  found  to  yield  satisfactory  results. 

In  order  to  compute  the  local  Mach  number  at  the  edge  of  the  boundary  layer 
at  *  =L,  a  linear  distribution  is  assumed  between  the  conical  shock  value 
and  the  value  at  x  =  l ,  v/hich  is  approximately  equal  to  2.0.  The  problem 
now  reduces  to  determining  where  conical  flow  is  reached. 

If  conical  flow  conditions  are  attained  on  the  body,  then  (Me ) i_  =  (Mc)Cone  • 

If  conical  flow  is  not  attained  on  the  body,  then  (Me)L  <  (Me)cone  in  accord¬ 
ance  with  the  assumed  linear  distribution. 

Employing  the  results  of  reference  25,  the  distance  required  to  attain 
conical  flow  on  a  blunt  cons  in  laminar  flow  is, 


•1 


It  should  be  noted  that  the  ”e"  conditions  are  evaluated  for  conical 
flow. 

Employing  the  subscript  "L"  to  designate  conditions  at  *  =  L,  and  the  sub¬ 
script  "c4' to  designate  conical  flow  conditions.  The  assumed  relation  for 
Me  at  x  =  L  is, 

L 

M_  -  2  +  (M.  -  2)  - 

L  c  ~c 

L  I. 

with  the  condition  that  if —  >  1,  set —  =  1 

c  c 

To  obtain  Te  at  x  =  L,  we  employ  the  relations, 


Solving  these  relations  simultaneously  and  assuming  that  bc  =  he 
we  obtain,  L 


where, 

hs  2 

1  ■»  1  +  0.2  for  y  *  1.4 

Tne  flow  properties  at  x  =  L  computed  in  the  above  manner  are  then  used 
in  the  previously  derived  equations  for  induced  pressure  drag,  pressure 
induced  skin-friction  drag,  and  transverse  curvature  induced  skin-friction 
drag. 

E.  Angle  of  Attack  Effects 

The  effect  of  angle  of  attack  is  to  increase  the  induced  drag  on  the  leeward 
side  and  decrease  Vne  induced  drag  on  the  windward  side.  Since  there  are 
no  rigorous  solutions  available  for  handling  this  problem,  the  approximate 
method  discussed  in  section  III  will  be  employed.  This  method  assumes 
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that  the  zero  angle  of  attach  solutions  are  applicable  when  the  flow  properties 
at  the  edge  of  the  boundary  layer  are  computed  using  tangent -cone  theory 
for  the  windward,  leeward,  and  side  meridians.  The  :,effective"  flow  prop¬ 
erties  are  then  obtained  by  averaging  the  properties  oi  the  four  meridians. 

“e  Pe  Te 

Therefore,  the  relations  given  in  section  III  for — ,  — ,  and— -will  be  util- 

U  p  T 

ized  in  the  induced  drag  equations.  “  00 

F.  Mass  Addition  Effects 


I  - 


The  effect  of  mass  addition  is  to  increase  all  the  induced  drag  contributions 
as  a  result  of  a  thickening  of  the  boundary  layer.  For  the  case  of  teflon 
ablation  on  an  eight  degree  cone  in  high  altitude  laminar  flow,  the  results 
of  reference  34  show  that  the  induced  pressures  due  to  viscous  interaction 
with  blowing  may  be  quite  significant. 

No  similar  calculations  have  been  made  during  this  drag  investigation,  since 
the  effects  of  mass  addition  were  not  to  be  considered  in  this  study.  How¬ 
ever,  the  procedures  followed  in  pursuing  induced  drag  calculations  with 
mass  addition  are  as  follows: 

a.  Employing  the  non-interacting  flow  properties  for  the  first  intera- 
tion,  compute  the  laminar  boundary  layer  displacement  thickness  deriva- 

tive  \-37-- )  fche  ratio  of  the  injected  mass  flow  to  external  mass 

flow  (  w  Vy  Y  This  is  obtained  using  laminar  similarity  theory  (Pro- 
\Pe  ue  / 

gram  1475)  {reference  29),  in  conjunction  with  the  previously  deter¬ 
mined  thermodynamic  and  transport  properties  of  the  injected  gas. 

b.  Employing  the  solution  derived  in  reference  35,  the  new  effective 
cone  angle  for  the  calculation  of  induced  external  flow  properties  is, 


c.  Employing  the  new  external  flow  properties  for  the  second  iteration, 
recompute  the  boundary  layer  parameters,  the  new  effective  cone  angle, 
and  the  corresponding  induced  external  flow  properties.  This  procedure 
converges  after  a  relatively  few  iterations. 

'The  final  requirement  would  be  to  correlate  the  results  of  the  solutions  into 
closed  form  expressions.  Since  this  tas.k  has  not  yet  been  undertaken,  it 
is  not  obvious  what  correlation  parameters  would  be  employed  for  this 
purpose. 


.  -a 
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V.  ERROR  ANALYSIS 


A.  Accuracy  Limits 

Sections  II,  III  and  IV  discuss  the  various  assumptions  and  appro:  imatious 
employed  in  the  development  of  the  equations  for  the  components  of  Cp.  This 
section  compares  the  values  of  total  Cp  predicted  by  these  equations  with  wind 
tunnel  and  ballistic  range  data.  In  addition,  a  comparison  of  predicted  Gp 
with  the  one  set  of  flight  data  available  (containing  accelerometer  and  pressure 
measurements)  is  made. 

The  accuracy  limits  for  total  Cp  predictions  are  30  percent  at  200,  000  feet 
and  10  percent  below  100,000  feet,  A  linear  interpolation  is  made  between  these 
two  altitudes.  As  will  be  shown  below,  these  accuracy  limits  have  been  met. 

B.  Data  Sources 

A  literature  search  was  conducted  at  three  libraries;  A vco  RAD,  Massachusetts 
Institute  of  Technology  and  thv  Redstone  Scientific  Information  Center.  A 
number  of  reports  were  evaluated  as  to  their  applicability  to  this  investigation. 
Those  not  used  are  listed  in  the  bibliography  with  a  note  indicating  the  reason. 
Those  used  are  references  36  through  41. 

C.  Typical  Drag  Coefficient  Data 

Typical  wind  tunnel  test  data  are  presented  along  with  ;he  predicted  values 
(solid  lines)  in  figures  24  through  27.  The  first  three  curves  show  Cp 
versus  a  for  several  Mach  numbers,  cone  angles  and  bluntness  ratios.  The 
predicted  curves  agree  well  with  the  data  except  at  the  larger  angles  of  attack 
where  they  tend  to  fall  below  the  measured  points.  This  deviation  is  greatest 
for  the  blunter  cones,  figure  26. 

Figure  27  shows  Cp  versus  Reynolds  number  for  a  cone  at  zero  angle  of 
attack.  AtR^*''  10^  there  appears  to  be  a  jump  in  the  experimental  values  of 
Cp  which  was  attributed  by  the  investigators  (reference  38  )  to  be  a  transition 
effect.  The  separate  laminar  and  turbulent  curves  shown  in  this  range  are 
drawn  through  the  values  predicted  for  each  ballistic  range  data  point.  These 
curves  indicate  that  some  of  these  test  conditions  might  have  been  partially 
turbulent. 

Figure  27  also  shows  a  trend  of  increasing  C  p  with  Mach  number  in  both  the 
predicted  curves  and  experimental  data.  This  trend  is  due  to  the  increase  in 
skin-friction  drag  and  induced  drag  as  the  Mach  number  is  increased.  To 
investigate  the  differences  between  theory  and  experiment,  the  ratio  of  pre¬ 
dicted  to  experimental  Cp  versus  Mach  number  is  plotted  in  figure  28. 

Here  it  is  seen  as  in  figure  27  that  the  predicted  values  are  somethat  less 
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Figure  28  DRAG  COEFFICIENT  RATIO  AS  A  FUNCTION  OF  MACH 


than  the  experimental  values.  No  curve  was  drawn  through  these  data  because 
of  the  scatter.  It  may  be  that  factors  other  than  Mach  number  play  a  part  in 
producing  the  variation  shown,  since  the  wall  temperature,  Reynolds  number, 
nose  bluntness  ratio  and  cone  angle  are  all  varying. 

Figure  29  shows  a  similar  ratio  versus  cone  angle.  Here  a  correction  trend 
is  indicated  but  no  line  is  drawn  due  to  the  number  of  parameters  involved  and 
the  data  scatter. 

D.  Error  Evaluation  Based  on  Ground  Tost  Data 

To  assess  the  reliability  of  the  methods  employed  for  evaluating  drag  coefficients, 
calculations  have  been  made  for  every  experimental  data  point  used  in  this 
investigation.  The  results  are  presented  in  table  H;  which*  includes  "the 
measured  total  drag,  the  ratio  of  predicted  to  experimental  drag  and  the  pre¬ 
dicted  values  for  each  drag  component. 

The  last  column  shows  the  minimum  effective  altitude  of  the  data.  It  is  required 
to  compute  this  factor  to  assess  whether  the  results  are  within  the  required 
accuracy  limits. 

With  a  simple  manipulation  the  expression 


derived.  The  bracketed  expression  (containing  density,  sound  speed,  and 
viscosity  of  air)  is  altitude  dependent  and  decreases  with  increasing  altitude. 
Fixing  the  left  side  of  the  equation  from  ground  test  parameters,  an  infinite 
number  of  combinations  of  altitude  and  length  will  satisfy  the  equation.  The 
larger  the  length  chosen  the  higher  the  effective  altitude.  It  is  clear  that  with 
a  length  range  of  1  to  12  feet  an  altitude  band  of  over  50,000  feet  could  be 
assigned  for  each  data  point.  To  assess  the  accuracy,  however,  the  smallest 
length,  1  foot,  was  chosen  to  establish  the  applicable  minimum  effective  altitude, 
hence  the  most  stringent  possible  accuracy  requirements. 

Figure  30.  presents  the  conversion  chart  used  to  establish  the  minimum  altitude 
for  each  data  point.  The  drag  coefficient  ratio  is  then  plotted  versus  this  altitude 
along  with  the  allowable  error  band  in  figure  31.  It  is  seen  that  close  to 
90  percent  of  the  points  fall  within  the  allowable  band.  It  is  considered  that 
this  verifies  the  accuracy  of  the  analytical  expressions.  There  does  appear 
to  be  a  trend  in  this  data,  as  in  the  previous  curves,  showing  that  the  pre¬ 
dicted  values  are  somewhat  less  than  the  experimental  values.  Additional 
analyses  are  required  to  determine  the  reasons  for  this  apparent  trend  and 
whether  adjustments  to  the  theory  are  warranted. 
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Figure  30  ALTITUDE  AS  A  FUNCTION  OF  REYNOLD5  NUMBER  FOR  VARIOUS  MACH 
NUMBERS  —  BODY  LENGTH  l  FOOT 
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:-E.  Applications  to  Flight  Conditions 

In  an  attempt  to  evaluate  the  accuracy  of  the  analytical  drag  equations  derived 
in  this  report,  the  equations  should  also  be  compared  to  flight  test  data. 

Applying  the  theory  to  flight  conditions  is  necessary  for  the  following  reasons: 

1)  The  primary  employment  of  the  equations  will  be  for  predicting  drag 
coefficients  in  reentry  flight. 

2)  Some  portions  of  the  analytical  expressions  were  based  von  flight 
environments;  u  e. ,  base  drag,  real  gas  curve-fits  to  the  reference 
enthalpy  method,  and  boundary  layer  transition. 

3)  The  ground  test  data  does  not  adequately  cover  the  Mach  number 
Reynolds  number  regime  encountered  in  this  study. 

4)  The  ground  test  data  does  not  consist  of  sufficient  turbulent  flow  data 
to  compare  with  theory;  also  it  is  not  known  how  much  of  the  boundary  layer 
is  turbulent. 

5)  Flight  test  data  is  presently  the  primary  means  of  evaluating  the 
effects  of  ablation  on  vehicle  drag. 

At  the  present  time,  there  is  a  very  small  amount  of  reliable  flight  test  data 
that  can  be  used  for  this  prupose.  The  reasons  for  this  scarcity  are  as  follows: 

1}  To  investigate  the  viscous  effects,  the  flight  vehicle  should  be  a 
slender  body  so  that  the  skin-friction  and  induced  drag  will  be  a  significant 
part  of  the  total  drag. 

2)  The  vehicle  should  be  instrumented  with  sensitive  accelerometers. 

The  results  of  integrating  the  accelerometer  measurements  should  yield  a 
trajectory  that  agrees  with  the  tracking  data. 

3)  It  is  also  desirable  to  have  a  vehicle  that  is  instrumented  with  pres¬ 
sure  taps  so  that  the  pressure  drag  portion  of  the  total  drag  may  be 
determined.  However,  on  slender  sharp  cones,  this  is  generally  not 
necessary  since  the  pressure  drag  at  high  Mach  numbers  can  be  calculated 
quite  accurately. 

At  the  present  time,  Avco  RAD  has  reliable  flight  test  drag  data  from  three 
conical  vehicles  which  are  very  nearly  pointed  {RN/RB  «  0.013).  The  data 
from  two  of  these  flights  is  presently  being  reduced;  the  results  of  one  of  the 
flights  have  already  been  documented  (reference  42),  The  total  drag  data 
from  this  flight  is  presented  in  figure  32,  together  with  the  results  of  the 
theory  which  has  been  developed  in  this  report.  The  estimate  of  the  altitude 
range  of  transition  was  obtained  from  the  results  presented  in  Section  VI. 
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Included  in  this  comparison  are  the  theoretical  results  obtained  when  the  effects 
if  mass  addition  are  neglected.  .The  actual  vehicle  was  covered  with  a  phenolic 
refrasil  heat  shield  (except  for  ‘Teflon  in  the  tip -region);  outgassing  began  at 
about  200,000  feet,  and  surface  ablation  commenced  at  about  100,  000  feet. 

The  skin-friction  drag  coefficient  was  reduced  for  blowing  by  the  approximate 
method,  employing  figure  lg' 

This  comparison  indicates  good  agreement  between  the  flight  test  data,  and  *he 
analytical  methods  presented  in  this  report;  the  differences  are  within  the 
limits  of  the  accuracy  requirements. 

Until  the  discrepancies  between  theory  and  test  data  can  be  more  thoroughly 
evaluated,  it  is  recommended  that  the  equations  derived  in  this  report  be 
employed  to  predict  the  drag  of  conical  flight  vehicles  to  an  accuracy  that  is 
within  the  specified  requirement.  One  of  the  reasons  for  the' differences 
between  theory  and  flight  test  data  might  be  the  method  used  to  account  for  mass 
addition  effects.  If  this  is  the  case,  than  it  would  be  necessary  to  employ  the 
mere  rigorous  methods  of  analysis  mentioned  in  sections  III- 5  and  IV-6. 


vi.  BOUNDARY  layer  transition  on  cones  ' 

The  objective  of  this  transition  study  was  to  develop  a  method  of  estimating  the 
altitude  at  which  a  conical  reentry  vehicle  (sharp  or  blunt)  would  begin  to  under¬ 
go  boundary  layer  transition.  In  order  to  show  how  such  a  method  has  been 
developed,  the  following  paragraphs  deal  with  a  discussion  of  the  transition  cri¬ 
terion,  how  the  criterion  i®  employed,  and  how  the  results  have  been  correlated. 

A.  Discussion  of  Transition  Criterion  - 

«■>—  ■  ■■  'I  . ...I  .1  1. 1.1  !■«.  -  II  

The  criterion  presently  employed  has  been  obtained  from  the  flight  test 
results  of  various  reentry  vehicle  programs.  Reference  43  provides  a  dis¬ 
cussion  of  the  flight  test  data,  how  it  was  interpreted,  and  how  it  wae  cor¬ 
related  to  yield  a  reasonable  criterion.  , 

This  investigation  of  flight  data  indicated  that  the  transition  Reynolds  num¬ 
ber  was  most  influe  need- by  the  local  Mach  number  at  the  edge  of  the  bound¬ 
ary  layer.  More  specifically,  correlation  was  achieved  when  the  computed 
values  of  local  Reynolds  number  at  the  observed  transition  altitudes  for 
each  flight  were  plotteci  versus  the.  corresponding  computed  values  of  local 
Mach  number.  These  results  are  represented  by  the  curve  in  figure  33 
which  is  a  "best  estimate"  of  the  band  of  data. 

The  data  actually  includes  largo  variations  of  the  other  parameters  which 
-  influence  transition:  nore  blunting,  angle  of  attack,  pressure  gradients, 
body  forcer,  Surface  roughness,  wall  temperature,  and  mass  injection 
rates-  Although  the  effect  of  each  of  these  parameters  is  known  to  be  sig¬ 
nificant,  the  net  effect  did  not  alter  the  observed  transition  altitudes  very 
~  much  (less;  than 35,  000  feet).  This  is  qualitatively  attributed  to  the  fact 
that  several  of  these  affects  tend  to  oppose  each  other.  For  example,  the 
low  wall  temperature  s  ye  suiting  from  a  lew  temperature  ablator  tends  to 
stabilize  the  boundary  tayer- however,  the  high  mass  injection  rates  associ¬ 
ated  with  low  temperature  ablators  tend  to  destabilize  the  boundary  layer 
(reference  44).  furthermore,  a  cool  wall  thins  the  boundary  layer,  magni¬ 
fying  the  destabilizing  effect  of  surface  roughness.  This  opposing  effect 
has  boon  Uxpeyimentaily  demonstrated.  (Reference  45  and  others). 

The  angle  of  attack  effect  is  generally  quite  significant,  and  it  is  an  effect 
which  is  moat  difficult  to  evaluate-  However,  based  on  the  flight  test  re¬ 
sults  shown  in  reference  43,  significant  differences  in  angle  of  attack  did 
not  alter  the  transition  altitude  by  more  than  10, 000  feet. 

The  effect  of  nose  blunting  is  partially  accounted  for  in  the  present  criterion 
since  the  actual  nose  bluntness  strongly  influences  the  local  Reynolds  num¬ 
ber  -  Mach  number. distribution  along  the  body.  The  consequence  of  nose 
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blunting  that  is  not  accounted  for  is  the  destabilizing  effect  of  the  adverse 
pressure  gradient  which  generally  follows  the  overexpansion  at  the  sphere  - 
cone  tangent  point.  For  low  bluntness  configurations,  this  effect  is  not 
very  critical  since  the  adverse  gradient  occurs  in  the  forward  region  of  the 
cone  where  the  Reynolds  number  is  low  compared  to  the  rear  portion. 

The  present  criterion  shown  in  figure  33  illustrates  how  significant  the 
local  Mach  number  effect  is,  covering  about  two  orders -of -magnitude  of 
(Rex)-j-g  or  about  100,  000  feet  in  altitude.  Although  this  graph  was  constructed 
from  flight  test  results  with  numerous  other  influencing  parameters  inherent 
in  the  data,  the  trend  is  quite  typical  of  most  of  the  wind  tunnel  data  where 
the  test  conditions  were  carefully  controlled  (reference  46  and  others). 

The  "bucket"  in  the  curve  at  Me  =  4  is  also  characteristic  of  wind  tunnel 
results. 

B.  Employment  of  Criterion 

In  employing  the  criterion  of  figure  33  to  predict  the  altitudes  where 
boundary  layer  transition  occurs,  it  is  necessary  to  evaluate  the  Reynolds 
number  -  Mach  number  (Rex  -  Me  )  distributions  along  the  edge  of  the  bound¬ 
ary  layer.  For  this  task,  the  following  methods  of  analysis  are  employed: 

a.  Sharp  Cone 

For  this  case,  M.  and  the  local  unit  Reynolds  number.  Re,  are  constant 
with  respect  to  surface  distance,  x.  The  values  of  these  quantities  as 
a  function  of  free -stream  conditions  and  cone  angle  are  then  easily 
obtained  from  conical  flow  theory:  i.e.,  reference  1. 

b.  Blunt  Cone 

For  this  case,  the  Rex  -  Me  distribution  along  the  edge  of  the  boundary 
layer  is  determined  from  the  blunt  body  viscid-inviscid  method  of 
analysis  (references  24>?/  Z5}  (Computer  Program  1115B).  This  pro¬ 
cedure  employs  method  of  characteristics  solutions  (reference  3) 
(Computer  Program  596)  for  the  pressure  distribution  and  shock  shape. 
The  blunt  body  viscid-inviscid  method  of  analysis  is  an  iterative  pro¬ 
cedure  as  outlined  below  for  a  particular  body  station: 

1)  Employing  the  known  pressure  (from  the  method  of  character¬ 
istics),  assume  a  value  of  entropy  for  the  external  flow.  These 
two  properties  are  then  used  in  Mollier  charts  to  obtain  all  other 
required  properties  external  to  the  boundary  layer. 

2)  Employing  these  external  properties  in  conjunction  with  the 
known  wall  conditions,  compute  the  mass  flow-  in  the  laminar  bound¬ 
ary  layer. 
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3)  By  performing  a  mass  balance  between  the  boundary  layer 
flow  and  the  free  stream  flow  which  is  entrained  in  the  boundary 
layer,  compute  the  stream -tube  radius  of  the  free -stream  mass 
flow, 

4)  Employing  the  known  shock  shape  (from  the  method  of  charac¬ 
teristics),  determine  the  shock  angle  at  the  location  where  the 
bounding  streamline  intersects  the  shock. 

5)  From  this  shock  angle  and  free-stream  Mach  number,  the 
entropy  behind  the  shock  is  determined.  This  is  compared  to  the 
assumed  entropy,  and  iterations  are  made  until  the  solution  con¬ 
verges. 

Some  typical  Rex  -  Me  distributions  resulting  from  the  blunt  cone  method  of 
analysis  are  presented  in  figures  34  to  36  for  cone  rngles  of  8,  iO,  and  15 
degrees  with  varying  lengths.  The  lines  of  constant  length,  x,  represent 
the  slant  length  of  the  equivalent  sharp  cone  (i.  e.,  surface  distance  from 
theoretical  sharp  tip  to  end  of  cone).  For  each  value  of  x,  the  correspond¬ 
ing  body  radius  is  given,  so  that  the  bluntness  ratio  may  easily  be  calcu¬ 
lated  for  a  specific  configuration.  The  actual  nose  radii  selected  corres¬ 
pond  to  a  bluntness  ratio  of  0.  05,  for  x  *  12  feet.  Additional  Rex  -  Me  dis¬ 
tributions  for  other  nose  radii  are  reported  in  reference  47. 

The  above  investigation  was  conducted  for  a  free-stream  velocity  of  25,  000 
ft/ sec  at  several  altitudes  encompassing  the  expected  range  of  transition. 

The  effect  of  decreasing  the  velocity  is  to  reduce  both  Rex  and  Me,  thus 
shifting  the  maps  down  and  to  the  left. 

The  distributions  presented  in  figures  34  through  36  show  how  Rex  and  Me 
increase  along  the  body  approaching  the  sharp  cone  value.  At  high  altitudes, 
conical  flow  is  reached  sooner  on  the  body  because  of  the  thick  boundary 
layer  which  is  capable  of  entraining  a  large  mass  flow.  At  lower  altitudes, 
the  body  length  required  to  reach  conical  flow  increases,  and  for  L  -  12  feet, 
the  results  shew  that  in  many  cases  conical  flow  will  not  be  attained.  The 
results  also  show  that  the  more  slender  the  cone  the  greater  the  distance 
required  to  reach  conical  flow.  A  closed  form  solution  for  the  distance  re¬ 
quired  to  attain  conical  flow  on  a  blunt  cone  has  been  derived  in  reference  25. 

In  using  these  Rex  -  Me  maps  to  determine  the  transition  altitude  and  the 
mode  of  transition,  the  criterion  shown  in  figure  33’  is  superimposed  on 
the  maps. 

This  procedure  yields  the  altitude  at  which  a  particular  body  station  (x)  will 
undergo  transition.  A  station  is  considered  to  become  turbulent  when  its 
R«x  -  Me  value  is  such  that  it  falls  just  above  the  criterion  curve.  By 
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determining  the  transition  altitude  for  sevjral  values  of  *,  the  mode  of 
transition  for  a  particular  vehicle  may  be  evaluated. 

"N 

C.  Correlation  of  Results 

The  transition  criterion  (figure  T3  and  the  (Rex-Me)  maps  (figures  34 
to  36,  plus  those  presented  in  reference  47)  have  been  employed  in  a 
generalised  study  in  order  to  correlate  the  transition  altitude  as  a  function 
of  cone  angle,  body  length,  bluntness  ratio,  and  free -stream  velocity. 

Based,  on  the  ranges  of  variables  considered,-  and  the  valid  ty  of  the  transi¬ 
tion  criterion  employed,  the  investigation  should  yield  reasonable  results 
for  the  following  conditions: 

cone  an^le:  8*  -  0  -  15* 

body  length:  lr  -  L  -  12' 

blmtness  ratio:  0  -  RN/RB  "  0r3 

free-stream  velo  '.ty:  18,000  -  -  25,000  ft/aec  (this  range  covers 

the  normal  cor'  dor  of  reentry  velocities)^ 

wall  temperature:  1,  000  -  Tw  -  5,  000* R 

angle  of  attack:  0*  -  a  -  6 

wall  material:  ablating  or  outgassing 

The  degree  of  accuracy  of  the  transition  altitudes  presented  is  expected  to 
be  witl  in  15,  0C0  feet  of  altitude. 

Non-ablating  vehicles  with  a  smooth  surface  finish  (i.e.,  beryllium)  are 
expected  to  undergo  transition  15, 000  to  20,000  feet  below  the  transition 
altitude  for  the  corresponding  ablating  /ehicle. 

The  generalized  study  was  performed  for  a  velocity  of  25,  000  ft;  sec  using 
the  methods  described  in  the  previous  section.  The  results  obtained  for 
specific  cases  of  body  length  and  nose  radiur  were  replotted  as  a  function 
of  bluntness  ratio.  These  plots  were  then  cross-plotted  at  bluntness  ratios 
of  0,  0.  1,  0.  2, "and  0.  3;  the  final  results  are  shown  in  figures  37  to  40. 

These  results  show  the  altitudes  where  transition  is  initially  expected  to 
occur  on  the  body.  Fc>r  all  the  configurations  investigated,  the  initial  sta¬ 
tion  to  undergo  transition  was  always  at  the  rear  of  the  body  (x  =L);  then 
as  the  altitude  decreases,  the  transition  location  proceeds  upstream.  In 
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the  case  of  low  bluntness  cona*  with  #<  8®,  the  R*x-  ^distribution  may 
show  that  transition  will  occur  first  in  the  forward  region  of  the  cone,  thus 
triggering  transition  ovet  the  major  portion  of  the  body. 

For  the  cases  shown  in  figures  37  to  40,  where  transition  occurs  first 
in  the  rear  portion,  there  are  two  types  of  information  that  can  be  obtained 
from  the  graphs:  1)  the  altitude  of  initial  transition  may  be  obtained  as  a 
function  of  6,  RN/Rg  ,  and  L. ,  and  2)  for  a  particular  vehicle  in  which  these 
three  parameters  are  known,  the  altitude  range  of  transition  may  be  deter¬ 
mined;  i.e.,  for  an  8-degree  cone,  12  feet  long  and  a  0.  10  bluntness  ratio 
at  a  velocity  of  25,  000  ft/sec,  figure  38  shows  that  transition  will  com¬ 
mence  at  130,000  feet  and  the  one  foot  station  will  become  turbulent  a S 
90,  000  feet. 

To  account  for  the  effects  of  free -stream  velocity  in  the  range  of  18,  000 
^  u„  -  25,000  ft/ sec,  the  results  shown  in  figures  37  to  40  may  be  used 
in  conjunction  with  a  velocity  correction  factor;  i.e. , 

Z-rjj  «  (Z-J-O )  +  AZ 

1  lR  uw  -  2^)00 

The  quantity  AZ  represents  the  increase  in  transition  altitude  as  the  velocity 
is  decreased  below  25,000  ft/ sec.  This  increase  in  altitude  is  most  signi¬ 
ficant  in  the  case  of  sharp  and  slightly  blunted  cones  having  high  values  of 
at  z  s  L.  The  reason  is  as  follows:  as  the  velocity  is  decreased.  Mc  and 
Rex  are  both  decreased.  However,  according  to  figure  33,  the  reduction 
in  SI,  causes  to  decrease  at  a  faster  rate  than  Rex.  Therefore, 

transition  occurs  earlier  when  the  velocity  is  decreased. 

Ia  the  case  of  highly  blunted  cones,  Me  ia  relatively  low  (2  <  84,  <  4)  and  as 
shown  in  figure  33  (R«x)xr  is  not  very  sensitive  to  Me  in  this  range. 
i.  refer®,  with  blunted  cones  the  transition  altitude  is  not  very  sensitive 
c  iree  stream  velocity. 

A  parametric  evaluation  cf  this  velocity  effect  is  shown  la  f  igure  41. 

These  results  represent  the  approximate  increase  iv  transition  altitude  as 
a  function  of  velocity  and  bluntness  ratio  for  an  average  cone  angle  and  body 
length. 

The  results  shown  .in  figures  37  to  41  represent*  a  first  attempt  to  cor¬ 
relate  transition  altitude  as  a  function  of  the  influencing  p*.ra.~neters.  The 
main  wea'ness  of  this  method  is  thit  it  fails  to  include  all  the  factors  which 
tend  to  destabilize  the  boundary  layer.  The  argument  presented  here  is 
that  the  local  Mach-numbe r  influence  on  the  transition  Reynolds  number  is 
the  dominating  effect,  while  the  net  effect  of  all  the  other  parameters  is 
relatively  small.  This  argument  stems  from  the  analysis  of  flight  test  data 
provided  in  reference  43, 


-94- 


VELOCITY  EFFECT  ON  TRANSITION  ALTITUDE 


To  check  the  accuracy  of  the  correlation  curves  presented,  they  have  been 
used  to  re-estimate  the  transition  altitudes  for  the  various  flight  tests  in 
which  reliable  transition  data  was  available.  In  all  cases,  agreement  was 
achieved  within  10,  000  feet  of  altitude.  In  predicting  the  transition  altitudes 
for  future  vehicles,  it  is  expected  that  these  results  will  be  valid  within 
15,000  feet  if  the  flight  conditions  are  within  the  specified  range  of  variables. 

As  more  test  data  becomes  available,  the  present  transition  criterion  will 
be  improved  to  include  more  of  the  influencing  parameters.  The  new  cri¬ 
terion  can  then  be  used  in  conjunction  with  the  {  Rex-Me)  maps  shown  in  this 
report  and  in  reference  47  in  order  to  develop  new  graphs  for  predicting 
the  transition  altitudes. 
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VIL  SPHERE  DRAG  AND  TRANSITION 


The  drag  coefficient  of  a  sphere,  like  that  of  a  cone,  is  the  turn  of  the  contribu¬ 
tions  of  pressure  drag,  base  drag,  skin  friction  drag,  and  induced  drag.  For 
spheres,  however,  the  pressure  and  skin  friction  drag  are  the  dominant  factors. 
The  induced  drag  for  spheres  is  negligible  for  all  cases  of  practical  interest. 

Base  drag,  although  not  negligible,  is  a  small  percentage  of  the  total  drag  at 
Mach  numbers  between  5  and  25. 

The  pressure  drag  is  usually  considered  to  be  a  function  of  Mach  number  only, 
while  the  skin  friction  drag  is  a  function  of  Reynolds  number.  Therefore,  the 
sphere  drag  coefficient  prediction  would  be  expected  to  have  the  form  Cp  «  f(M) 

+  s(R«L>.  Data  which  confirm  this  expression  from  references  48  through  54 
are  presented.  No  data  on  wall  temperature  effects  were  found. 

A  plot  of  drag  coefficient  versus  Mach  number  is  shown  in  figure  42.  These 
data  correspond  to  flight  below  150,  000  feet,  based  on  the  sphere  sire  and  flight 
conditions  called  for  in  this  study.  The  predicted  constant  value  of  Cp  «  0.  92 
fits  the  data  well  within  a  ±  10  percent  band  for  M>  3.  It  is  assumed  that  cn  re¬ 
mains  constant  through  M=25  although  the  highest  Mach  number  tested  was  approxi 
mately  16.  Test  results  from  other  blunt  bodies  substantiate  this  assumption. 

Drag  coefficient  versus  Reynolds  number  is  plotted  in  figure  43.  Values  in 
the  region  10*  £.  Re  £  10^  represent  Mach  numbers  as  low  as  2.5.  Here,  too, 
the  constant  value  of  C3  =  0.  92  fits  the  test  data  within  a  ±  10  percent  band  for  a 
simulated  flight  altitude  of  almost  200,  000  feet,  discribed  by  a  3.  75-inch  diam¬ 
eter  sphere  at  Mach  5.  At  K„L5xlG*  or  an  altitude  greater  than  150,000  feet, 

CD  increases  with  decreasing  Rc.  A  curve  fit  of  this  trend  yields 

500 

C0  «  0.92  ♦  r - 

S. 

This  equation  is  applicable  throughout  the  range  of  interest. 

B,  Boundary  Layer  Transition 


For  spheric?!  configurations,  a  limited  parametric  study  has  been  made  to  cor¬ 
relate  transition  altitude  as  a  function  of  sphere  radius, R,  ,  and  free-stream 
velocity. 

For  this  investigation,  a  constant  value  of  transition  Reynolds  number  {r*x)jr  - 
300,  000  was  assumed.  To  determine  the  local  Reynolds  number  distribution 
along  the  sphere,  the  fluid  properties  were  computed  using  a  Newtonian  pres¬ 
sure  distribution  and  the  entropy  behind  a  normal  shock.  The  maximum  unit 
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SPHEXF.  D". AG  COEFFICIENT  AS  A  FUNCTION  OF  SEYNOEOS  NUM.ESt 


Reynolds  number  occurs  at  the  sonic  point  where  the  mass  flow  is  maximum. 
Employing  the  computed  unit  Reynolds  number  distribution  in  conjunction  with 
the  surface  distance  from  the  stagnation  point  ansi  the  assumed  value  of  (&*,)!*• 
the  transition  altitudes  have  been  determined.  The  results  are  presented  in 
figure  44  for  velocities  o i  5,000,  15,  COO,  and  25,  000  ft/sec.  Thes*  graphs 
show  that  the  great*;  the  velocity,  the  higher  the  transition  altitude,  Thif 
occurs  because  the  local  Reynolds  number  increases  with  velocity  wliile  tht 
local  Mach  number  is  essentially  constant. 
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VIn  CONCLUSIONS 

The  analytical  expressions  presented  in  this  report  satisfy 
the  basic  requirements  of  the  study: 

1)  All  equations  are  in  terms  of  free  stream  properties, 
body  geometry,  and  surface  material. 

2)  The  accuracy  requirements  have  been  met  based  on 
comparison  with  extensive  wind-tunnel  an (1  ballistic 
range  data,  and  one  set  of  flight  test  data. 

3)  The  specified  ranges  of  variables  have  been  covered. 

The  additional  requirement  that  has  been  satisfied  is  the 
determination  of  the  transition  altitude  ran^e  as  a  function 
of  the  several  influencing  parameters 

The  equations  that  are  presented  in  this  report  for  cones  have 
been  programmed  for  digital  computer  solut;on  (Program  1954) 
(Reference  55).  The  results  of  this  program  have  been  compared 
to  more  rigorous  solutions  and  there  has  been  excellent  agreement. 
The  program  is,  however,  limited  to  the  ranges  of  variables 
covered  in  this  report. 

In  order  to  improve  and  expand  the  drag  solutions  presented 
herein,  it  is  suggested  that  the  following  additional  tasks  be 
perfoi  med: 

1)  Employ  more  rigorous  solutions  in  order  to  determine  the 
effects  of  blowing  on  skin  friction  drag  and  induced  drag. 

2)  Extend  the  Mach  number  and  cone  angle  range  lo  include 
the  transonic  region  and  cone  angles  as  shallow  at  4°. 

3)  Develop  drag  equations  for  the  triconic  shape  as  the  general 
configuration,  noting  that  cones,  bi -conics,  and  cone-cylinder- 
flares  are  special  cases  of  the  triconic. 

These  three  tasks  can  be  performed  employing  existing  capabilities 
The  additional  tasks  that  require  highly  advanced  state-of-the-art 
techniques  would  be  involved  with  the  determination  of  all  the  drag 
components  at  angles  of  attack  greater  than  the  cone  half-angle.  The 
approximations  employed  in  this  investigation  for  viscous  effects  at 
o( j/o  should  be  valid  only  at  small  angles  of  attack. 
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